Spin Polarization and Transport of Surface States in the Topological Insulators Bi2Se3 

and Bi2Te3 from First Principles 



Oleg V. Yazyev, Joel E. Moore, and Steven G. Louie 

Department of Physics, University of California, Berkeley, California 94720, USA and 
Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 

(Dated: December 7, 2010) 

We investigate the band dispersion and tiie spin texture of topologically protected surface states 
in tlie bulk topological insulators Bi2Se3 and Bi2Te3 by first-principles methods. Strong spin-orbit 
entanglement in these materials reduces the spin-polarization of the surface states to ~50% in 
both cases; this reduction is absent in simple models but of important implications to essentially 
any spintronic application. We propose a way of controlling the magnitude of spin polarization 
associated with a charge current in thin films of topological insulators by means of an external 
electric field. The proposed dual-gate device configuration provides new possibilities for electrical 
control of spin. 
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The recently discovered three-dimensional topological 
insulators (TIs) realize an unconventional electronic 
phase driven by strong spin-orbit interaction (SOI). The 
striking feature of these heavy-element materials is the 
existence of metallic Dirac fermion surface states char- 
acterized by an intrinsic spin helicity: the wavevector of 
the electron determines its spin state. A net spin den- 
sity is thus produced upon driving a charge current at 
the surface of a TI. Considerable effort has been devoted 
recently to possible applications of this property of TIs 
to novel electronic devices, e.g., to spintronics 01 and 
topological quantum computing 

Understanding the properties of bulk TIs has mostly 
relied on few-band phenomenological models which treat 
the helical surface states as fully spin-polarized. How- 
ever, one has to keep in mind that SOI entangles the 
spin and orbital momentum degrees of freedom thus re- 
ducing spin polarization. This is especially so for the cur- 
rently investigated bulk TIs in which SOI is of electron- 
volt magnitude due to the presence of bismuth, the ele- 
ment with the largest atomic number Z which has stable 
isotopes d, Q- 

In this Letter, we apply first-principles methodology to 
study the topological surface states in Bi2Se3 and Bi2Te3, 
the "second generation" materials which are currently 
considered as reference TIs (1, Q . This approach is free of 
empirical parameters and treats topological surface states 
on an equal footing with other states across the electronic 
spectrum. The degree of spin polarization of the topo- 
logical surface states is found to be significantly reduced 
in both bismuth materials, and this reduction will affect 
most spintronic applications of these materials. We use 
our results to propose a way of controlling the degree 
of spin polarization associated with the charge current 
in thin TI slabs by applying a transverse electric field. 
The proposed spintronic device is compared to a related 
proposal based on spin-polarized edge states in graphene 
[lO| and to conventional semiconductor quantum wells. 



First-principles electronic structure calculations have 
been performed within the density functional theory 
(DFT) framework employing the generalized gradi- 
ent approximation (GGA) to the exchange-correlation 
functional 11]. Spin-orbit effects were treated self- 
consistently u sing fully relativistic norm-conserving pseu- 
dopotentials [1^ acting on valence electron wavefunc- 
tions represented in the two-component spinor form [l^ . 
A plane-wave kinetic energy cutoff of 35 Ry has been 
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FIG. 1: (Color online), (a) Band structure of the 5 QL slabs 
of Bi2Se3 and Bi2Te3 (lines) superimposed with the projected 
band structure of the corresponding bulk materials (shaded 
areas), (b) Evolution of the surface state band dispersion in 
the vicinity of V point as a function of slab thickness. Shaded 
areas show the bulk band gaps, (c) Band gaps at the F point 
induced by the interaction between the surface states as a 
function of slab thickness. 
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employed for the wavefunctions. Surface states in Bi2Se3 
and Bi2Te3 bulk TIs were investigated using (111) slab 
models of varying thickness (1-5 quintuple layers (QL) 
equivalent to '^1-5 nm) with atomic positions taken from 
experimental data using a supercell geometry Ii4L Such 
an approach has already been used previously [1, [l^ [3 ■ 
The effects of applied electric field were simulated by 
adding a sawtooth-like potential profile with a constant 
slope in the slab region [17|. The PWSCF code of the 



QUANTUM-ESPRESSO distribution [Till was used in the 
present study. 

We start our discussion by considering the band struc- 
ture of the bulk materials and slabs of finite thickness. 
Figure [TJa) shows the band structure of 5 QL slabs su- 
perimposed with the bulk band structure projected onto 
the surface Brillouin zone (BZ). The band structures of 
the slab and bulk systems are aligned by matching the 
potential in "bulk-like" region of the slab to the bulk po- 
tential. Within the DFT-GGA approach the band gaps 
of Bi2Se3 and Bi2Te3 are 0.31 eV and 0.08 eV, respec- 
tively. The value for Bi2Se3 is in good agreement with 
another theoretical report (0.32 eV [l9[) and experimen- 
tal data (0.35 eV [20|), a fortuitous result considering 
the well-known tendency of DFT to underestimate band 
gaps [2l'|. The calculated band gap for Bi2Te3 is consis- 
tent with another first-principles result ^22jand is smaller 
than 0.165 eV measured experimentally [23|. 

The topological surface states are distinguished in the 
slab band structures as the "Dirac cone" feature at the 
r point, especially clear in the case of Bi2Se3. The Dirac 
point energies (below, — eV is set for conve- 
nience) are situated 0.04 eV and 0.10 eV below the corre- 
sponding bulk valence band maxima (VBM) within DFT- 
GGA. In slab geometry, the surface state bands are dou- 
bly degenerate with corresponding surface states of the 
opposite spin helicities localized at the opposite surfaces. 
The topological surface bands in bulk TIs do not respect 
electron-hole symmetry, and the Fermi velocity vp shows 
rather complex behavior. At the VBM energy the calcu- 
lated vp = 4.6 X 10^ m/s in Bi2Se3 agrees with experi- 
mental results Q . The magnitude of v-p increases with in- 
creasing energy and reaches 6.4xl0^m/satthe bulk con- 
duction band minimum (CBM) energy for wavevectors 
along the F-K direction. However, along the F-M direc- 
tion there is a noticeable softening {v-p — 3.8 x 10^ m/s) 
due to the hexagonal warping effect. The surface-state 
band of Bi2Te3 shows stronger anisotropy in agreement 
with experimental observations of a strong hexagonal 
warping At E ^ 0.1 eV, vp = 4.1 x 10^ m/s 

(3.3 X 10^ m/s) along the F-K (F-M) direction. The 
variation of vp along the F-K direction is weak, while 
along the F-M direction it reduces to 1.2 x 10^ m/s at 
the GBM energy. 

Figure[T]Jb) illustrates the evolution of the surface band 
dispersion as a function of slab thickness. Vacuum poten- 
tial alignment is used for the comparison of the slab sys- 
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FIG. 2: (Color online), (a-c) Expectation values of the spin 
operators {Sa{k)) {a — x,y,z) for the topological surface 
states along the selected path (A — F — B) in momentum space 
(d) for 4 QL slabs of 612863 and Bi2Te3. Electron (solid lines) 
and hole (dashed lines) surface state bands are distinguished. 
The surface normal corresponds to z direction. Coordinate 
system is right-handed (d, inset). 
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In thin slabs, the interactions between the states 



localized at the opposite surfaces opens a gap {Eg ) at the 
Dirac point. The magnitude of the gap decays rapidly 
with increasing slab thickness [Fig. [Ijc)], and the values 
are in quantitative agreement with other first-principles 
results na. The dispersion of the surface state bands is 
essentially converged by 3 QL and 4 QL in Bi2Se3 and 
Bi2Te3, respectively. On the other hand, in thin slabs, 
quantum confinement pulls apart the bulk valence and 
conduction band edges thus increasing the energy range 
in which topological surface states are separated from 
bulk-like states. 

The striking feature of the surface states of TIs is 
the helicity of the spin polarization vector P{k) — 
{2/h)[{Sx{k)), {Sy{k)), {Sz{k))] along the constant energy 
contours. Here, the expectation values of spin oper- 
ators {Sat{k)) = {h/2){i(;i{k)\aa\ipi{k)) {a = x,y,z), 
where ipi{k) are the two-component spinor wavefunc- 
tions, ctq the corresponding Pauli matrices. The re- 
cent spin- and angle-resolved photoemission spectroscopy 
(spin-ARPES) experiments have indeed indicated a one- 
to-one locking of the momentum and the direction of spin 
polarization vector pointing along {k x z) 1241 . First- 
principles results have confirmed this picture [25|. How- 
ever, we would like to stress that the magnitude of spin- 
polarization P{k) can be reduced from the maximum 
value of 100% since the electron spin quantum number is 
no longer conserved in systems with SOI [26]. Spin-orbit 
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coupling leads to Bloch states that are not separately 
spin eigenstates but instead have entanglement between 
their spin and orbital parts. An illustration of the ef- 
fects of SOI on the spin projections in the case of bcc 
iron is given in Ref. |27l Such effects are expected to 
be especially pronounced in bismuth materials where the 
exceptionally strong SOI, of the order of an electron-volt 
0, 13 ' would make the surface states a mixture of bulk 
states from a broad energy range. 

Figures [2ja)-[2Uc) show the calculated expectation val- 
ues of the spin operators {Sa{k)) for both electron and 
hole states of the surface band in Bi2Se3 and Bi2Te3 along 
X and y directions in momentum space [see Fig. [2][d) for 
definition]. The surface normal corresponds to z direc- 
tion and the standard convention of the right-handed co- 
ordinate system applies [Fig. [^Jd), inset]. We find that 
the spin polarization vectors are aligned preferentially in 
the xy plane and the helicity is left (right) handed for the 
conduction or above (valence or below E^) bands. 
This is consistent with previous experimental measure- 
ments [2^ and theoretical calculations 25]. However, 
the magnitudes of the spin projections are significantly 
reduced with respect to the nominal value of ±fi./2. Close 
to the r point we find ~50-60% spin-polarization for 
both materials. Its magnitude tends to decrease with in- 
creasing energy. The magnitude of spin projections close 
to the Dirac point appears to be rather insensitive to the 
slab thickness (not shown here). We have studied the mi- 
croscopic origin of the reduction of spin polarization and 
find that the surface Bloch state is a complex superposi- 
tion of different spin directions on different atoms; hence 
the reduction of the spin polarization cannot simply be 
understood in a tight-binding picture built from atomic 
eigenstates of total angular momentum. A remarkable 
out-of-plane spin projection {Sz{k)) develops in the F- 
K direction as it was recently ascribed to the hexagonal 
warping effect (28| . In general, we find somewhat larger 
reduction of the spin polarization and more pronounced 
hexagonal warping effects on the out-of-plane spin pro- 
jection in the case of heavier Bi2Te3. 

From the point of view of technological applications, 
an attractive feature of the TI materials is the intrinsic 
spin polarization of the current carried by the topolog- 
ical surface states. For charge carriers (electrons) mov- 
ing along y direction (ky > 0) the transferred spin po- 
larization (i.e. spin polarization of the current injected 
through a transparent tunneling barrier) is along x axis 
and its geometric average {Px{E)) = (7r/4)|P(_E)|. We 

J dkVy{k)P^{k)S(E-Eg} 



define {P:,{E)) 



where Vy{k) is 



J dkvyik)S{E-Eg) 

the carrier velocity along the transport direction. In the 
discussed materials {Px{E)) w 50% and shows only weak 
dependence on E, but changes its sign when the char- 
acter of charge carriers (electrons vs. holes) changes. 
Moreover, in TI films the opposite spin helicities of the 
charge carriers at the opposite surfaces would result in 
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FIG. 3: (Color online), (a) Illustration of the effect of electric 
field Eaxt applied across a thin slab of bulk TI. The filling 
of the two surface state bands and the corresponding spin 
textures for electrons and holes are shown, (b) Splitting of 
the surface state bands by applied electric field {5E — 0.1 eV, 
solid lines) in a SQL slab of Bi2Se3 compared to the band 
dispersion at zero field (dotted lines). The inset shows a close- 
up view of band gaps at the charge neutrality point, (c) 
Individual (surfaces 1 and 2) and total ballistic conductances 
in y direction per device width, Gy, and (d) spin polarization 
along X axis, Px, for SE = 0.1 eV and SE = eV as a function 
of chemical potential fi {fi = eV corresponds to the charge 
neutrality point). The symbols have same meaning in (c) and 
(d). 



zero net spin polarization. 

We anticipate that improved control over spin trans- 
port can be achieved in few-nanometer thick slabs of 
Bi2Se3 and BioTeB which can be produced by molecular- 
beam epitaxy [29l . Isoj , vapor deposition [3l|, H^l , or exfo- 
liation [3 31 . Tuning the chemical potential Mjn such thin 
slabs by gating has already been achieved [321 ] . Similarly 
to bilayer graphene [s^, the TI slabs are characterized 
by imperfect screening due to the semimetallic nature of 
the surface states. Thus, in the dual-gate device configu- 
ration or by coupling TI thin films to a ferroelectric sub- 
strate it is also possible to control splitting 5E — E^—Ef 
between the Dirac point energies E^ and E^ of the two 
surface bands given with respect to the common chemical 
potential [Fig. [3l^a)]. Such splitting of ~0.1 eV produced 
by the SiC substrate has been reported for Bi2Se3 thin 
films [2^. Close to the Dirac point, the response of a 
thin TI film of thickness d to an applied electric field 
Ecxt is expected to be dominated by the bulk screening 
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contribution, 5E = \Ecxt\d/£\\{0). The static dielectric 
constant along c-axis e||(0) = 75 for Bi2Te3; probably 
smaller for Bi2Se3 [s^. Figure EJb) shows the surface 
state band dispersion calculated from first-principles for 
a 3 QL slab of Bi2Se3 under applied electric field Ecxt re- 
sulting in 6E = 0.1 eV (here, Ecxt = 1-72 V/nm for the 
supercell dimension along z axis equal to 4.36 nm). The 
conical surface state bands shift rigidly with respect to 
each other while the bulk-like states (not shown) remain 
practically unaffected by the applied field. Figure [3jc) 
shows the individual surface state Sharvin conductances 
along y direction (i.e. the number of ballistic channels 
per device width along x axis given in the units of lattice 

1 2 

constant a), Gy and Gy, and their sum Gy. 

The applied field lifts the semi-metallic conductance 
close to the charge neutrality point. Note that for both 
6E = 0.0 eV and SE = 0.1 eV there is a 0.006 eV band 
gap opening at the charge-neutrality point due to the 
hybridization of the surface states localized on the op- 
posite sides of the thin film. For E^ < fj. < the 
charge carriers are electrons (holes) at surface 1 (2). Im- 
portantly, in this range of fi the transferred spin polar- 
ization along X direction is of the same sign for both 
surface conduction channels [Fig. [3l[d)]. The total trans- 
ferred spin polarization (Px) is non-zero in the case of 
applied electric field. We note that these conditions are 
away from the symmetric situation that has been argued 
to lead to exciton condensation [36]. The magnitude of 
(Px) can be tuned by changing fi and achieves maximum 
its value {Px) = 0.608 at ^ = The sign of {Px) 

can be changed by simply inverting the direction of E^xt ■ 
That is, the dual-gate device based on a thin slab of TI 
permits to control independently both the conductance 
and the degree transferred spin polarization of injected 
current in a broad range, certainly a highly attractive 
opportunity for spintronics applications. 

Such spintronic device can be compared to another 
theoretical proposal which realizes electric control over 
spin transport. This device based on zigzag graphene 
nanoribbons relies on the peculiar spin-polarized edge 
states of topologically-trivial origin [lOj. While in TI 
slabs the degree of transmitted spin-polarization is lim- 
ited by the SOI and by the geometric average prefactor 
7r/4, the localized states at the zigzag edges of graphene 
are always 100% spin-polarized. On the other hand, the 
degree of spin-polarization can be continuously tuned in 
the TI slabs. The direction of P is fixed with respect to 
the current direction at the TI surfaces, while graphene- 
based magnetic systems are characterized by low mag- 
netic anisotropy as a result of very weak SOI in carbon 



current also induces a spin density. The basic difference 
is that the Rashba coupling leads to two Fermi surfaces 
with opposite spin directions at each momentum, and 
the resulting current-induced spin densities nearly can- 
cel if the SOI is small. A single TI surface has only one 
spin state at each momentum and hence no cancellation. 
Quantitatively, we write the Rashba SOI as ah{k x a)z 
where a has units of velocity and a are Pauli matrices. 
The cancellation between Fermi surfaces reduces the spin 
density by a factor of order a/vp compared to a model TI 
surface with a single Dirac cone, where wp is the Fermi ve- 
locity. Taking one example, a in a InGaAs/InAlAs quan- 
tum wells was measured to be 7-14 x lO'^ m/s [1^, more 
than an order of magnitude less than up ~ 5 x 10^ m/s 
in Bi2Se3 and Bi2Te3 [l,!!!. 
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and by relatively short spin correlation length [37, 

In a more general context, it is instructive to com- 
pare the spin transport effects at the TI surfaces to those 
in other 2D systems, e.g. to the corresponding effect 
induced by Rashba SOI in a 2D electron gas, where a 
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